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ABSTRACT 

Globular clusters are source of 7-ray radiation. At GeV energies, their emission is 
attributed to magnetospheric activity of millisecond pulsars residing in the clusters. 
Inverse Compton scattering (ICS) of ambient photon fields on relativistic particles dif- 
fusing through cluster environment is thought to be the source of GeV to TeV emission 
of globular clusters. Using pair starved polar cap model 7-ray emission from synthetic 
millisecond pulsar was modelled. In addition to pulsar emission characteristics, the 
synthetic pulsar model yielded spectra of electrons escaping pulsar magnetosphere. 
To simulate 7-ray emission of globular cluster, both products of synthetic millisecond 
pulsar modelling were used. Gamma-ray spectra of synthetic millisecond pulsars resid- 
ing in the cluster were summed to produce the magnetospheric component of cluster 
emission. Electrons ejected by these pulsars were injected into synthetic globular clus- 
ter environment. Their diffusion and intera ction, both, with clus t er ma gnetic field and 
ambient photon fields, were performed with lBednarek k, Sitarekl (|2007t ) model yielding 
ICS component of cluster emission. The sum of the magnetospheric and ICS compo- 
nents gives the synthetic 7-ray spectrum of globular cluster. The synthetic cluster 
spectrum stretches from GeV to TeV energies. Detailed modelling was preformed for 
two globular clusters: Tcrzan 5 and 47 Tucanae. Simulations are able to reproduce 
(within errors) the shape and the flux level of the GeV part of the spectrum observed 
for both clusters with the Fermi/LAT instrument. The synthetic flux level obtained in 
the TeV part of the clusters' spectrum is in agreement with a H.E.S.S. upper limit de- 
termined for 47 Tuc, and with emission level recently detected for Ter 5 with H.E.S.S. 
telescope. The synthetic globular cluster model, however, is not able to reproduce the 
exact shape of the TeV spectrum observed for Ter 5. 

Key words: gamma-rays: theory - globular clusters: individual (Ter 5, 47 Tuc) - 
pulsars: general - radiation mechanisms: non-thermal 



1 INTRODUCTION 

Recently, globular clusters (GCs) were established as a 
source of GeV emission. So far, the LAT instrument onboard 
the Fermi Space Telescope has detected 15 globular clusters 
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(|Abdo et al.ll2009al . l20ld iKong et aLllJOlol : iTam et aLllJOllI ) 
emitting photons in the energy range from ~ 200 MeV to 

~ 50 GeV. All of them are seen as po int sources. 

For 8 of the 7-ray emitting GCs JAbdo et al.ll2010t ) the 
placement of the emission is coincident with the position 
of the clusters. In these cases the observed 7-ray emission 
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does not vary in time. Moreover, their differential photon 
spectra can be fitted with a power law with an exponential 
cut-off The derived spectral indices are F < 2.0, and the 
cut-off energies £c fall in the range between 1 GeV and 4.5 
GeV. However, for two of these GCs (M62 and NGC 6652) 
the d etermined spectral cut-offs are not statistically signif- 
icant (jAbdo et al.ll201Cl ). Nonetheless, the observed charac- 
teristics of the 7-ray emission of globular clusters resembles 
that observed for millisecond pulsars (MSPs) d etected with 
the Fe rmi/LAT instrument at GeV energies (JAbdo et al.l 
[2009b|). Thus, the population of MSPs residing in the clus- 
ters can be responsible for their 7-ray radiation. 

iTam et al.l (|201ll ) reported detection of 7 more 7-ray 
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emitting globular clusters. Interestingly, their 7-ray emission 
is slightly shifted with respect to the position of the cluster. 
The extreme case here is NGC 6441 for which the maximum 
of the GeV e mission falls far o utside the cluster tidal radius 
(see fig. 4 of iTam et al.ll201lh . The 7-ray emission of these 
7 GCs does not vary i n time. How e ver, i n contrast with 
the clusters detected bv I Abdo et al.l (|201(J ). their observed 
photon spectra are best fitted by a power law. The exception 
here is the spectrum of NGC 6441 tha t is better describ ed by 
the exponentially cut-off power law (jTam et al.l 120111 1. For 
two of the clusters (Liller 1 and NGC 6624) the 7-ray spectra 
extend to energies above 40 GeV. The power-law character 
of these spectra and their energy extent cannot be solely 
attributed to the magnetospheric emission of MSPs residing 
in these clusters. Possibly, the observed emission from these 
two GCs results from the inverse Compton scattering (ICS) 
of the ambient photon fields on relativistic electrons injected 
into globular cluster environment by MSPs. 

In addition to the Fermi/LAT detection in the GeV do- 
main, the 7-ray emission at the TeV energies was detecte d 
towards Terzan 5 globular cluster (|Abramowski et al.ll2011f ). 
It is slightly extended beyond the GC tidal radius. More- 
over, the peak of the emission is signifi cantly shifted with 
respe ct to the cluster centre (see fig. 1 of lAbramowski et al.l 
I2OIII ). The observed spectrum can be approximated by a 
power law with a photon spectral index F ~ 2.5. The ICS 
photons may substantially contribute to the observed TeV 
radiation from Terzan 5. However, the ICS origin of TeV 
emis sion cannot account for the observed source morphol- 
ogy. lAbramowski et al.l (|201l|) give different possible expla- 
nations of the TeV radiation detected from the direction of 
Ter 5. Source coincidence is one of the explored possibilities. 
In such case the TeV 7-rays could originate from a pulsar 
wind nebula associated with a radio-quiet pulsar. However, 
the probability of such scenario is very low. Other possi- 
bility is that the TeV emission could have hadronic origin 
(for details see lAbramowski et al.ll201ll ). This scenario could 
explain the power-law emission spectrum. No TeV emission 
from the dire ction of other globular clus ters has been de- 
tected so far jAnderhub et a l. 2009; Aha ronian et al.|[2009l : 
lMcCutcheonll2009l). 

iHui et al.l (|201l|) studied these 15 globular clusters de- 
tected in the 7-rays JAbdo et all l2009al . I2OIOI : iKong et al] 
I2OI0I : iTam et al.ll201lh in search of correlations between the 
cluster 7-ray luminosity L-y and various cluster properties. 
Such correlations could possibly facilitate explanation of the 
origin of the 7-ray emission of these GCs. iHui et al.l (|201lh 
investigated the correlation between L^ of the GC and two- 
body encounter rate Fc, metallicity [Fe/H], absolute visual 
magnitude, and photon densities of the infrared and the op- 
tical Galactic background. The weakest correlation is found 
between L-y and the absolute visual magnitude of the clus- 
ter. However, strong correlation is found between L^ and 
the two-body encounter rate and also the metallicity of the 
cluster. The reported co rrelations are co nsist ent with re- 
sults p reviously found bv lHui et al.l (|2010|) and lAbdo et al.l 
(|201CI ). Moreover, they find that the cluster 7-ray luminosity 
increases with the increase of the density of the background 
photon fields, both the optical and the infrared ones. The 
correlation of L^ and Fc suggests that the 7-ray emission 
of globular clusters is related to the population of objects 
formed through dynamical interactions in the cluster. This 



in turn leads to the conclusion that L^ depends on the num- 
ber of MSPs residing in the cluster JAbdo et al.ll2010l ). 

Metal-rich globular clusters are more likely to contain 
bright low-mas s X-ray binaries (LMXBs). As proposed by 
llvanoval (|2006r ). this can be related to the different stel- 
lar structure of main sequence donors with masses between 
~ O.85M0 and ~ 1.25M0. Metal-poor main-sequence stars 
do not have an outer convective zone, while metal-rich ones 
do. Absence of the convective zone shuts down magnetic 
braking, which in turn makes orbital shrinkage in these bi- 
naries less efficient. Thus, preventing them from becoming 
bright X-ray sources, while opposite is the case for binaries 
with metal-rich main-sequence donors. As low-mass X-ray 
binaries are thought to be progenitors of MSPs, higher prob- 
ability of LMXBs formation for metal-rich globular clusters 
translates into a higher formation rate of MSPs for these 
clusters. Thus, the observed correlation between L~, and 
[Fe/H] can link the origin of the 7-ray emission in GCs to 
the population of MSPs. 

The observed correlation between L^ and the density of 
the photon fields cannot be accounted for if the 7-ray radia- 
tion would originate solely from the magnetospheric activity 
of MSPs. However, such correlation is expected when the 7- 
rays are produced via ICS. Then the emitted power in ICS 
process is directly proportional to the energy den sity of the 
ambient photon field. The corre l ations found by Hui et al.l 
(i20Tl] ) (and also lHui et al.ll2010l : JAbdo et aLllJOlOl ) suggest 
that the interplay between the number of MSPs and the 
density of the background photon fields determines the 7- 
ray luminosity of globular clusters. 



2 THEORETICAL APPROACH TO 7-RAY 
EMISSION OF GLOBULAR CLUSTERS 

In the past years, 7-ray emission from globular clusters was 
modelled assuming it results from either combined magneto- 
spheric radiation of millisecond pulsars residing in a cluster, 
or from IC scattering of different photon fields on relativistic 
electrons diffusing through GC interior. 

The latter scenario was studied in detail by 
iBednarek fc Sitarekl (|2007l . hereafter BS07). In their model, 
BS07 assumed that relativistic electrons originate either in 
magnetospheres of pulsars populating a globular cluster or 
are accelerated at the shock waves arising from collisions 
of pulsar winds. In the first case BS07 assumed monoen- 
ergetic spectra for the electrons; in the second case they 
used power-law spectra. In both cases, after being injected 
into the cluster environment the relativistic electrons dif- 
fuse through the cluster interacting with its magnetic field. 
Moreover, on their way out electrons up-scatter photons a- 
rising from stellar population in the cluster and photons of 
the cosmic microwave background (CMB). The modelling of 
BS07 predicts 7-ray spectra of GCs spanning from GeV to 
TeV energies. 

Similarly to BS07, ICheng et all (|2010l ') modelled 7-ray 
emission from globular clusters as arising from ICS processes 
taking place in the globular cluster environment. However, 
it was assumed that the relativistic electrons originate only 
from magnetospheres of millisecond pulsars residing in the 
cluster. The injected electron spectra were monoenergetic . 
The relativistic leptons in the model of ICheng et al.l (|2010|) 
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up-scatter not only the stellar photons and the CMB, but 
also soft photons originating in the Galaxy (i.e. infrared 
background and stellar photons from the galactic disk). The 
resultant 7-ray emission spans GeV to sub-Te V energies. It 
is imp ortant to point out that the model of ICheng et al.l 
(J201CI 1 predicts rather high sub-TeV flux for globular clus- 
ters. Moreover, it also predicts that the 7-ray emission 
should extend beyond the cores of GCs (a radius > 10 pc). 



IVenter et al.l (J2009f ) assumed that both processes, ICS 
and magnetospheric 7-ray emission of millisecond pul- 
sars, contribute to the overall 7-ray brightness of globu- 
lar clusters. Under assumption that an accelerating elec- 
tric field in the magneto s pheres of MSPs is unscre ened 
(JMuslimov fc Tsvganill99^ : lHarding fc Muslimovll 19981 ) due 
to insufficient pair creati on, 7-ray emis s ion fr om MSPs was 
simulated. In their work, IVenter et al.l ((20091) used popula- 
tion of cluster MSPs as a sole source of relativistic electrons 
in the GC environment. The energy spectra of electrons in- 
jected into the cluster were computed using the same model 
of the pulsar magnetosphere as the one used for simulating 
7-ray emission of MSPs. Using Monte Carlo method, av- 
erage cumulative magnetospheric spectrum and cumulative 
electron injection spectrum originating from 100 MSPs were 
calculated. This allowed for calculating GeV-to-TeV spectra 
for globular clusters 47 T uc and Ter 5. Simil arly to BS07, 
to produce ICS component I Venter et al.l (|2009l ) up-scattered 
the stellar photons and the CMB photons. 



We note that while modelling of GeV-TeV spectral char- 
acteristics of GCs puts constraints on the possible number 
of MSPs populating the cluster and on the cluster mag- 
netic field (given energy den sities of ambient ph oton fields 
and diffusion coefficients; e.g. IVenter et al.ll2009l ). X-ray en- 
ergy domain yields additiona l constraints on physica l prop- 
erties of GC. In particular, iBuesching et al.l (|2012l ) mod- 
elled synchrotron emission maps at radio and X-ray energies 
for Terzan 5. Comparison of the modelled X-ray emission 
map with the de tected diffuse X-ray emission from Ter 5 
(JEger et al.ll2010l l allowed to put additional constraints on 
diffusion process of relativistic particles within GC. For more 
details on the content of GCs, their high energy obser vations 
and modelling see recent review bv lBednarekl (J201ll ). 



In this paper we present the results of the numerical 
modelling of 7-ray radiation from globular clu s ters. In our 
study it is assumed (similarly as in lVenter et al.l (|2009f )l that 
the cluster 7-ray emission results from combined magneto- 
spheric activity of millisecond pulsars residing in the clus- 
ter, and from ICS scattering of ambient photon fields on 
relativistic electrons propagating through the cluster. These 
electrons are injected into the cluster environment by the 
MSPs. Detailed description of the synthetic millisecond pul- 
sar model used for computing 7-ray spectra of pulsars resid- 
ing in the cluster and also spectra of electrons ejected from 
their magnetospheres can be found in Sect.[3l In Sect.[4]de- 
tails of calculations of synthetic 7-ray spectra of globular 
clusters are presented. Comparison of the simulations with 
the existing observational data for selected globular clusters 
can be found in Sect. (5] Conclusions are given in Sect.[6l 



3 NUMERICAL MODEL OF SYNTHETIC 
MILLISECOND PULSAR 

In this section we discuss models of 7-ray production 
in the open magnetosphere (distances comparable to the 
light cylinder radius) of pulsars, a nd motivate why the 
pair s tarved polar cap model (PSPC; iMuslimov fc Hardind 
l2004bl ) is preferred by us for the MSPs residing in globu- 
lar clusters. We also discuss results of PSPC model that are 
used to calculate magnetospheric and ICS components to 
the 7-ray GeV-TeV emission of synthetic globular cluster 
(Sect. Ell. 



3.1 Justification for the selected synthetic 
millisecond pulsar model 

In the standard polar c ap model approach (see e.g., 
iHarding fc Muslimovlll998l ) the accelerating electric field is 
being screened out at distances close to the neutron star 
surface (~ few stellar radii), which means that the acceler- 
ation gap fills only a small fraction of the open magneto- 
sphere. Such description may be applicable to classical pul- 
sars (with ages Tc < 10^ yr) where in their magnetospheres 
pair production through interaction of curvature and/or in- 
verse Compton scattered photons with the magnetic field is 
sufficient to co mpletely screen the electric fi eld at low alti- 
tudes (see e.g.. iMuslimov fc Harding l2004a ). However, the 
situation is completely different for much older pulsars, in- 
cluding millisecond ones. For these pulsars the primary elec- 
trons may keep accelerating and at the same time keep emit- 
ting high energy photons up to very high altitudes without 
significant pair production (see e.g., Muslimov fc Hardind 
l2004bl : lHarding et al.|[2005[ ). lHarding et al.l (|2005D show that 
for millisecond pulsars the primaries can be accelerated up 
to distances comparable to the light cylinder radius. 

The characteristics of some of the 7-ray light curves 
of the millisecond pulsars and their behaviour with respect 
to the radio ligh t curves (i.e. misal ignment of 7-ray peak 
and radio peak; lAbdo et al.l l2009bl ) seem to favour outer 
magnetosphere as a region where the 7-rays originate. How- 
ev er, the 7-ray and r adio light curve modellin g performed 
bv I Venter et al.l (|2009l ') and I Venter et al.l (|2012l ') showed that 
the group of 7-ray bright MSPs is non-uniform in terms of 
the preferred emission model. Three subclasses can be dis- 
tinguished - two of the subclasses favour the outer magne- 
tosphere as the place of origin of the 7-ray radiation (either 
in the framework of the outer gap or the slot gap model). 
However, there are MSPs whose 7-ray emission can only be 
explained invoking the exten ded polar cap model, i.e. the 
pair starved polar cap model (|Muslimov fc Hardind l2004br ) 
where the 7-ray emission is produced in the whole volume 
of the pulsar open magnet osphere. 

In the slot gap model (JMuslimov fc Hardindl2004al ) the 
approximate solution to the accelerating electric field E^ | suf- 
fers from a sign reversal. In the out er gap model non-dipolar 
magnetic field (|Cheng et al.l I2OO0I ) needs to be invoked in 
order to obtain copious electron-positron pair creation. In 
turn, the created pairs are needed to control the gap dimen- 
sions. The pair starved polar cap model, similarly to the slot 
gap model, suffers from the sign reversal in the solution for 
the unscreened accelerating electric field i?|| (see Sect. 3.3). 
Nonetheless in the past years it was widely studied (see e.g.. 
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IVenter et al.ll2009l . and references therein) in the context of 
7-ray emission of millisecond pulsars. 

Because properties of MSPs in globula r clusters are dif- 
feren t from the ones of the field population (|Camilo fc Rasid 
[2OO5I) it is justified to suspect that the 7-ray emission mecha- 
nism operating in the magnetospheres of these globular clus- 
ter MSPs may be different from the one operating in the clas- 
sical pulsars and/or in the field millisecond pulsars. Thus, 
in the presented study the pair starved polar cap (PSPC) 
model is chosen to simulate the 7-ray emission of millisec- 
ond pulsars. 



3.2 Basic assumptions of the synthetic minisecond 
pulsar model 

To simulate the 7-ray emission of millisecond pulsars a 3 
dimension al (3D) num erical model of pulsar magnetosphere 
was used (|Dvkg 120021 ). The model operates within a space 
charge limited flow framework. The magnetic fi eld of a pul- 
sar is in the form of a retarded vacuum dipole (jDyks et al.l 
[2OOJ) for which the curvature radii of magnetic field lines 
are deter mined in the inertial f rame of reference. Moreover, 
following iDvks fc RudakI (|2002l ') the special relativity effects 
like aberration and time-of-fiight delays are treated accord- 
ingly throughout the calculations. 

Important change introduced into the numerical code 
was an inclusion of a more realistic description of the elec- 
tric field -B|| (for details see Sect.|33J. It is assumed that 
acceleration of particles takes place in the whole volume de- 
termined by the last open magnetic field lines. The elec- 
trons are injected at the polar cap surface at the Goldreich- 
Julian rate, and they have low initial energy. As there is no 
screening of the electric field, like in the sta ndard polar cap 
model (see e.g.. lHarding fc Muslimovll 19981 ). acceleration of 
the particles takes place even at distances comparable to 
the light cylinder radius. Such model is called a pair starved 
polar cap model (|Muslimov fc Hardinejl2004bl ). The imple- 
mented electric field takes into account the general relativis- 
tic effect of dragging of inertial fram es (|Muslimov fc TsvganI 
ll992l : lMuslimov fc Hard"in3l2004bt ). 



3.3 Structure of the accelerating electric field 

In the 3D numerical model of the millisecond pulsar magne- 
tosphere the accelerating electr ic field E \\ in the form pro- 
posed by iMuslimov fc Hardind (|2004a ) was implemented. 
The field itself extends up to the light cylinder and consists 
of three different prescriptions (_Ei, E2 and i^s) determining 
B|| at different heights above the pulsar surface. 

In the region close to the neutron star surface (further 
referred to as the near regime) the acc elerating electric field 
is desc ribed by the formula derived bv lMuslimov fc Tsygad 
(|l992h under assumption that the electric field must disap- 
pear at the stellar surface (^n = 0): 
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system (r, 9, 0), where r is the radial distance from the 
neutron star centre, 6 is the magnetic colatitude, and <j} is 
the magnetic azimuthal coordinate. This is also the case for 
all the _B|[ formulae presented in this Section. The angle a 
in Eq. [T]is the angle between the rotation and the magnetic 
axis, 5 = 6/6{ri) is the magnetic colatitude scaled with the 
polar angle of the last open magnetic field line, 77 = r/Rt^s is 
the radial distance scaled with the star radius, and z = rj—1 
is the altitude above the neutron star surface scaled with the 
star radius. The electric field is derived using a small-angle 
approximation, in which the polar angle of the last open 
magnetic field line ^(77) <^ 1. This results in: 
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Here Oq is the magnetic colatitude (polar angle) of the base 
of the line at the stellar surface. Other quantities used in 
Eq. [T]are: 

f^-RNS _ el _ 2GMns 
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where Mns is the mass of a neutron star, e = rg/Rj^s is 
the compactness parameter (r^ = 2GMns/c^ is the gravi- 
tational radius of a neutron star), and k is the magnitude 
of the general relativistic effect of the frame dragging at the 
stellar surface measured in the stellar angular velocity fl. 
The function Jm is a Bessel function of order m. The ki and 
ki are the positive roots of the Bessel functions Jo and Ji . 
Functions 7^ and "/i are defined as; 
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Formulae for f(ri), HJrj) and 5{r]) were first defined in 
iMuslimov fc Tsvgad (|l992h and can be found therein. The 
near regime solution Ei is applicable for altitudes rj—1 <C 1. 
At moderate heights (6q -^ rj ~ 1 <C -^ — ; further 
referred to as the moderate regime ) the e lectric field is given 
by eq. 14 of lHarding fc Muslimo-JI (|l998l 'l: 
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Close to the light cylinder (further referred to as the 
far regime) the elec t ric fie ld is described by eq. 35 of 
IMuslimov fc Hardind (|2004bl ): 
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where A — iif(l)/-\//(l), and 77c i s a radial parameter 
at wh ich the electric field saturates (|Muslimov fc Hardind 
l2004al ). 

In order to obtain a smooth transition between the elec- 
tric field working in the moderate regime (Eq.[5| and the one 
working in the far r egime (Eq. [S]) the formula proposed by 
IMuslimov fc Hardin g (2004b) is used: 



The above equation is expressed in the magnetic coordinate 
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(7) 
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The value of 77c para meter is determined for each magnetic 
field line separately. iMuslimov fc Harding (|2004b|) showed 
that for millisecond pulsars the best-fit value of rjc falls in 
the range ~ 3 — 4. Moreover, they point out that for a given 
pulsaj spin period the parameter rjc is a function of an in- 
clination a and an azimuthal coordinate at the polar cap. 
In the numerical code before electrons are injected into 
the pulsar magnetosphere, we find the accelerating electric 
field solution in the whole volume of the open magnetosphere 
up to the light cylinder distances. When finding the full 
-B|| solution that could accelerate particles also in the dis- 
tant parts of the magnetosphere (so a combina tion of Equa- 
tionsQ][S]and[n]with the use of Eq.O, following I Venter et al.l 
l|2009f ). we require that the resultant electric field: 

- should be negative for all ry in the range between 1 and 
rjLC, where 77LC = r/Rhc\ 

- should transit from the near to the moderate E\ | regime 
(matching Ei with E2) when Ei ~ E2; 

- should transit from the moderate to the far E\ \ regime 
(matching E2 with £3) smoothly, which is achieved by ap- 
plying Eq. [7] and appropriate matching of ric ; 

- should converge to Es, for the distances close to the light 
cylinder. 

For the implemented electric field (Eq. [l] [5] and (Gjl a sign 
reversal of the electric field occurs at some altitude. This 
happens especially for those magnetic field lines that cross 
the null-charge surfacqj within the volume of the open mag- 
netosphere limited by the light cylinder radius. When E\\ 
changes sign, the particle oscillations occur. In order to alle- 
viate this somewhat problematic situation it is required that 
JS|I field is negative for all 77. If this condition is not fulfilled 
for a certain magnetic field line, the line is left out from the 
calculations so no acceleration of particles takes place along 
this line. To obtain smooth transition from E2 to E^ we 
take a br oader range of possib l e value s of r]c than the one 
found bv iMuslimov fc Harding! (|2004bf ) for millisecond pul- 
sars. The values of r]c are chosen from the range between 1 
and 6. In the matching procedure the largest possible value 
of r]c that ensures the smooth transition between the electric 
field solutions given by Eqs. [S]and[S]is chosen. 

The distribution of the rjc parameter across the polar 
cap for the case of Bs = 10* and different values of P and a 
is presented in Fig. [T] When the inclination increases a re- 
gion with progressively smaller values of rj^ increases at the 
polar cap side facing away from the rotation axis (<^ = 180°). 
This is true for all presented spin periods. For a — 80° large 
part of the magnetosphere is excluded from the calculations 
(no solution of _E|| is found). Such behaviour results from 
the fact that when we increase pulsar inclination the mag- 
netic field lines situated in the part of the magnetosphere 
at large distances from the rotation axis start crossing the 
null-charge surface. Thus, in order to mitigate the problem 
of particle oscillations the open magnetic field lines crossing 
the null-charge surface are excluded from the calculations. 

Previously, the 7-ray emission of millise cond pulsars 
was m odelled by e.g. IVenter fc De Jagen (|2005l ) and lZaiczvkl 
(|2008l ). Their calculations were carried out in the framework 



^ null-charge surface is the surface in the pulsar magnetosphere 
separating volumes filled with the particles of opposite charge. 



a - 20° «!>„ = 90° 



a - 20° -p^ = 90° 





-1.0 -0.5 0.0 0.5 1.0 



a - 50° ■t^ = 90° 



-1.0 -0.5 0.0 0.5 1.0 



a - 50° «i>^ = 90° 




11 -1.0 -0.5 0.0 0.5 1.0 'H -1.0 -0.5 0.0 0.5 1.0 





-1.0 -0.5 0.0 0.5 1.0 




-1.0 -0.5 0.0 0.5 1.0 



Figure 1. The distribution of ric parameter across the pulsar 
polar cap. The images show results for synthetic MSP with Bs = 
10* G. The left column shows the case with P = 1.5 ms, while 
the right column shows the case with P = 9.5 ms. The inclination 
changes along columns: top row a = 20°, middle row a = 50°, 
and bottom row a = 80° . Different values of the parameter rjc are 
colour-coded (see colorbars for the exact values). White regions 
are those for which no -Eii solution is found. Small and big solid 
circle mark the distance across the polar cap equal to 0.5rpc and 
l.Orpc, respectively. Here rpc is the polar cap radius. Horizontal 
red dashed line shows = 0° while vertical one points to = 90° 
{<j> is the magnetic azimuth angle). 



of the extended polar cap model. However, they considered 
accelerating electric field i5|| of the form given by Eqs.[T]and 
[5] Thus, electrons were acce lerated only up to the distances 
^ T^LC- IVenter et al.l ((20091) performed the calculations of 
the 7-ray characteristics of the millisecond pulsars using the 
pair starved polar cap model, thus including also Eq.[6]in the 
description of the accelerating electric field. The conditions 
that need to be met when finding the solution of E\ | along 
the magnetic field lines are similar for IVenter et al.l (|2009 ) 
and the implementation of the PSPC model presented in 
this work. However, one difference can be found. In order 
to match the el ectric field solutions in the near - moderate 
distance regime IVenter et al.l 1)20091 ) used the approximate 
formula 77b ~ 1 + 0.0123P~°'''''', where rjh is the distance at 
which the transition between Ei and E2 takes place. In this 
case, the distance rjh is independent of the position at the 
polar cap, so the change from Ei to E2 for all open magnetic 
filed lines takes place at the same distance from the neutron 
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star surface. In the presented study, however, the distance 
of transition between E\ and E2 is determined for each of 
the magnetic field Unes independently. The transition takes 
place at the distance where the condition of E\ ~ E2 is ful- 
filled. Thus, the transition height from the near to moderate 
regime varies across the polar cap. This in turn results in dif- 
ferent transition heights from the moderate to far regime, 
so the sli ghtly different r?c v alues, with respect to the ones 
found by I Venter et al.l (|2009|). The general behaviour of the 
?7c parameter with the increasing inclination is similar in 
both cases. Additionally, the step in the numerical grid of 
the Tjc values used in the -E|| calculations has direct impact 
on the exact distribution of 77c across the polar cap. Thus, the 
final result is the o utcome of the two e ffects. The observed 
differences between IVenter et al.l (|2009l ) and the presented 
implementation of the PSPC model, though observed, have 
minor effect on the calculated 7-ray emission characteristic 
of the synthetic millisecond pulsars. 

Having found the accelerating electric field structure 
across the pulsar magnetosphere, electrons are injected at 
the polar cap with low initial energy (70 = 10). The par- 
ticles are distributed across the polar cap [j] — 1) with 
the Goldreich- Julian charge de nsity that includes general 
relativistic corrections (see e.g.. lMuslimov fc Hardindl 19971 : 
iHarding fc Muslimov|[l99l ): 
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where a = \J\ — e/ri is the red-shift function. Their start- 
ing positions coincide with the foot-points of the magnetic 
field lines. Then the primary particles are being followed in 
their motion along the magnetic field lines. Moving through 
the pulsar magnetosphere the particles loose their energy 
through curvature radiation being at the same time acceler- 
ated by the electric field. 

The emitted curvature photons interact with magnetic 
field. However, secondary pairs (electrons and positrons) are 
produced in insufficient numbers to screen out the electric 
field. The majority of curvature photons escape pulsar mag- 
netosphere without being absorbed. At the same time, infor- 
mation on their energy and directio n of motion is collected to 
produce photon maps (|Dvkal2002f ). which give information 
on pulsar emission characteristics (the number of photons 
dN^ radiated per unit time, per unit energy d£, into unit 
solid angle dfi). The photon maps are used to obtain pulsar 
spectra and light curves for different observing angles C,. 

Numerical calculations were performed for discrete 
distribution of pulsar parameters (P, Bs, a). Pul- 
sar spin period was selected from the set of P = 
(1.5, 2.0, 3.5, 5.0, 6.5, 8.0, 9.5) ms. Magnetic field values were 
chosen from the set of Bs = (1.0, 3.5, 4.5, 5.5, 7.5, 9.0) x 10* 
G. Pulsar inclination angle was chosen from a range between 
10° and 80° with a step of 10° . Such parameter selection al- 
lowed for obtaining 336 different synthetic millisecond pul- 
sar cases. For each modelled case a 7-ray photon mapo and 



spectrum of electrons ejected from pulsar magnetosphere 
(Sect. 13. 4p were obtained. 



3.4 Spectra of electrons ejected from millisecond 
pulsar magnetosphere 

In the calculations we also collect information on the Lorentz 
factors 7ei of the electrons escaping the millisecond pulsar 
magnetosphere. Examples of the electron ejection spectra in 
the form logj^g(d/V'ci/d_Bci) versus logj^Q 7ci are presented in 
Figure (2] The dA/'ei is the number of primary electrons of 
given energy Ef.\ = 7ei rrieic^ (rriei is the mass of the elec- 
tron) escaping pulsar magnetosphere per unit time. For the 
case of Bs = 10* G and a — 20° the electron ejection spectra 
are single-peaked (the left column of Fig. \^. The Lorentz 
factors of the particles contributing to the peak are in the 
range ~ 10*'^ — lO''. For these cases, with the increase of 
the spin period from 1.5 ms to 9.5 ms the spectral peak 
becomes flatter and shifts towards the lower values of 7ei. 
For P = 1.5 ms the peak centre is at 7ei ~ lO®*^, while for 
the case of P = 9.5 ms it is at 7ei ~ 10^". With the in- 
crease of the inclination a pronounced low energy tail starts 
to develop, which eventually for slowly rotating synthetic 
millisecond pulsars (e.g. P = 5 ms and 9.5 ms) evolves to 
form a separate low energy peak. Such double-peaked elec- 
tron ejection spectrum can be seen for the case of P = 9.5 ms 
and a = 50° (the middle plot in the bottom row of Fig. (2]). 
Similarly as for the cases with a — 20° , with the increasing 
P the shift of the whole spectrum towards the lower ener- 
gies is observed for the cases with a — 50° and 80°. For 
the model with P — 9.5 ms and a = 80° this results in the 
fact that the spectrum spans the range of 7ci between ~ 10* 
and ~ 10^. The similar behaviour in the electron spectra is 
also observ ed for the mo dels with larger values of the mag- 
netic field (|Zaiczvkll2012l ). However, the development of the 
low energy tail followed by the emergence of the low energy 
peak takes place much later in terms of the pulsar inclina- 
tion (higher a) and spin period (larger P) than for the low 
Bs pulsars. 

The overall behaviour of the electron ejection spectra is 
linked to the electric field structure -E|| in the pulsar magne- 
tosphere (see Sect. 13.3]) . For a given case with P and Bs the 
development of the low energy peak in the electron spec- 
trum starts when more i5|| solutions characterised by the 
low values of the 77c parameter are found across the mag- 
netosphere. This in turn results in the fact that primary 
electrons accelerated in different parts of the pulsar mag- 
netosphere reach progressively lower energies at the light 
cylinder. When the electric field structure (determined by 
the value of r]c) evolves into the state when within approx- 
imately half of the magnetosphere high rjc electric field so- 
lutions are found, and the other half is characterised by the 
very low values of r/c or no solutions at all, then the double- 
peaked characteristics of the electron ejection spectrum sets 
in. Such situation happens exactly when the inclination be- 
comes > 40° (Fig. [1]). For a given case with Bs and a as 
the pulsar period increases the magnitude of the electric 



^ photon map yields information on pulsar emission characteris- 
tics (the number of photons radiated per unit time, per unit en- 
ergy, into unit solid angle) as a function of pulsar rotation phase 



ip and viewing angle (^ llDvksll2002l) . Pulsar emission spectrum for 
selected f is obtained through integrating strip of the photon map 
A^ wide over the pulsar rotation phase. 
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of the radiation observed for GCs (|Abdo et al.ll2009al . I2OI0I : 
iKong et al.ll2010l : iTam et al.ll201ll ) seem to confimi this sce- 
nario. The relativistic electrons injected into the cluster en- 
vironment by the MSPs can up scatter the ambient photon 
fields (cosmic microwave background, infrared background, 
and starlight either originating from stellar population in the 
globular cluster or from stars in the Galaxy). This inverse 
Compton scattering scenario has been proposed to explain 
the GCs 7-ray emission in the GeV and TeV energy range 
llBednarek fc SitareklboOTi : ICheng et al.|[201ol : I Venter etall 
uOOgf ) . The observational properties of the 7-ray emission of 
GCs cannot be explained solely with the magnetospheric ac- 
tivity scenario or the ICS scenario. Thus, most probably the 
interplay between the two models gives rise to the overall 
GeV and TeV emission of the globular clusters (see Sect.[T] 
for more details). 




Figure 2. Spectra of electrons escaping through light cylinder for 
the model with _Ba = 10* G. The F-axis in each of the subfigures 
shows logio(A/V'ei/rf-Eel [s'-'MeV-l]). Spin period P of 1.5, 5.0 
and 9.5 ms is presented in the top, middle and the bottom row, 
respectively. Inclination of the pulsar a changes from 20° , 50° to 
80° from the left column to the right one. 



field decreases. Thus, on average the primary electrons reach 
smaller energies at the light cylinder and the whole ejection 
spectrum moves towards lower values of 7ei. Finally, when 
the pulsar magnetic field strength is increased, the average 
magnitude of the electric field in the pulsar magnetosphere 
is increased. The previously described double-peaked spec- 
tral behaviour is also present for high Bs cases. However, it 
becomes pronounced only for slowly rotating (P > 5 ms) 
and highly inclined {a > 50°) synthetic pulsars because of 
the much larger average magnitude of the electric field. For 
the same reason, the spectrum shift towards lower electron 
energies is far less severe for the high magnetic field cases 
[Bb = 9 X 10* G), where the low energy end of the spectrum 
shifts only to 701 ~ 10"". For the low magnetic field cases 
[Bb = 10* G) the low energy end of the ejection spectrum 
can shift even to 7ei ~ 10"*. 



4 GAMMA-RAY EMISSION OF SYNTHETIC 
GLOBULAR CLUSTER 

Globular clusters have recently been established as a source 
of the 7-ray radiation (see Sect.[l|. Their high energy emis- 
sion is attributed to the ensemble of millisecond pulsars re- 
siding in their cores. Combined magnetospheric activity of 
the millisecond pulsar population is thought to be the source 
of the cluster emission above 100 MeV. The spectral shapes 



4.1 Numerical scheme of the synthetic globular 
cluster 

In order to simulate the 7-ray emission of a synthetic glob- 
ular cluster both the magnetospheric activity of the popu- 
lation of millisecond pulsars residing in the cluster, and also 
the ICS scenario is taken into account. The magnetospheric 
component (Sect.|4]2} is constructed using the results of the 
pair starved polar cap model calculations for different cases 
of millisecond pulsars (Sect. I3.3|l . The ICS component is 
calculated in two steps, first a cumulative spectrum of elec- 
trons injected by ensemble of synthetic millisecond pulsars 
into the cluster environment is computed. The cumulative 
electron spectrum is constructed using the electron ejection 
spectra calculated with the pair starved polar cap model (see 
Sects. 13. 4|) . Later, the cumulative spectrum is fed into the 
numerical model of the cluster (BS07) to produce the ICS 
radiation (Sect.|33J- 

The numerical procedure that calculates the 7-ray spec- 
trum originating from the ensemble of millisecond pulsars 
residing in the core of the synthetic cluster, and also the 
cumulative spectrum of electrons ejected from the magne- 
tosphcrcs of those synthetic MSPs is written in the IDL 
language. 

It is assumed that the number of millisecond pulsars re- 
siding in the synthetic cluster is A^psr. These pulsars reside 
in the cluster core, which in the numerical calculations is re- 
solved assuming that these pulsars are concentrated in one 
point placed in the cluster centre. In the first step, the basic 
parameters (a spin period P, a magnetic field strength Bs 
and an inclination angle a) are selected for each of the mil- 
lisecond pulsars residing in the modelled cluster. The values 
of these parameters are randomly selected from the discrete 
distribution determined by the synthetic MSP population in 
the database (Sect. lTS)) . The ra ndom selection is per formed 
with the IDL function randomu (|Park fc MilleiJIlQSSl '). 



4.2 Magnetospheric emission of the population of 
millisecond pulsars in the cluster 

Having selected basic parameters (P, Bs, a) that charac- 
terise each of A'psr millisecond pulsars populating our syn- 
thetic globular cluster, we are ready to compute the magne- 
tospheric contribution of the ensemble of MSPs to the high 
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energy emission of the synthetic cluster. Each of the selected 
MSPs emits high energy photons originating in the radiation 
processes that take place in the pulsar magnetospheres. As 
already discussed in Sect. O in the PSPC model only cur- 
vature photons are emitted. Within the framework of 
the studied model and the selected range of pulsar 
parameters (P, Bb, a) there is no additional con- 
tribution (with respect to the curvature radiation) 
at the GeV energy range from syn chrotron radia- 
tion d ue to secondary pairs (see e.e:.. lRudak &: Dvksi 
I1999I ). Moreover, any synchrotron radiation due to 
cyclotron resonant absorption of radio waves is be- 
yond the scope of the studied model. 

To create the cumulative magnetospheric contribution 
from the millisecond pulsars, first a viewing angle C, had to be 
selected for each pulsar. Its value is randomly selected from 
the sin C, distribution. Having selected the viewing angle, the 
observed 7-ray spectrum is constructed for each millisecond 
pulsar populating the modelled cluster: 



10"° r 
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where (^ is the randomly selected viewing angle, tp is the 
pulsar rotation phase and _D is a distance to the simulated 
cluster. This in turn allows for computing the overall mil- 
lisecond pulsar magnetospheric contribution to the observed 
7-ray emission of the synthetic cluster: 



AfpSR 



(10) 



In the left column of Fig.|3]examples of the cluster 7-ray 
spectra resulting only from the magnetospheric emission of 
population of millisecond pulsars are presented. In the top 
left panel of Fig.|3]one thousand cluster 7-ray spectra in the 
form of f -FWpgii(^) ''■'"^ shown. It is assumed that 20 millisec- 
ond pulsars reside in the modelled cluster. The distance to 
the cluster is 9.6 kpc. The obtained spectral characteristics 
of the synthetic cluster is diverse both in terms of the photon 
cut-off energy and the maximum flux density level. The en- 
ergy of the cut-off f cut falls in the range between 10^ MeV 
and ~ 5 X 10* MeV. The cluster spectra characterised by 
the lowest value of £cut are also the spectra with the low- 
est maximum flux density level ~ lO"'^'^ erg~^ s~^ cm~'^. 
The highest generated maximum flux density level is ~ 
10"^" erg~^ s~^ cm"'^. Assuming that there are only 20 
pulsars in the cluster we are able to generate the magneto- 
spheric 7-ray spectra of the globular cluster that differ by 
two orders of magnitude in terms of the flux density level. 
In the bottom left panel of Fig. [3] three hundred cluster 7- 
ray spectra are displayed. The simulated globular cluster 
harbours A'^psa — 100 millisecond pulsars, and is at the dis- 
tance of 9.6 kpc. The cut-off energies of the simulated mag- 
netospheric cluster spectra span the energy range between 
~ 2 X 10^ MeV and ~ 2 X 10* MeV. The maximum level of the 
flux density varies between 



3 X 10"" erg"^ s"^ cm~^ and 



3 X 10~^° erg~^ s"'^ cm~^. The comparison of the models 
presented in the top left {NpsR ~ 20) and in the bottom 
left (ApsR = 100) panel of Fig. |3] shows that the synthetic 
globular clusters harbouring smaller number of millisecond 
pulsars can produce similar magnetospheric 7-ray spectra 
(in terms of the cut-off energy and the flux density level) as 
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Figure 3. The magnetospheric 7-ray spectra (left column) and 
the cumulative electron spectra [right column) computed for dif- 
ferent realisations of synthetic globular cluster. The magneto- 
spheric 7-ray spectra are presented in £F^pgj^{£) form (i.e. a 
7-ray flux described by Eq. 1101 multiplied by the photon energy 
£). The cumulative electron spectra are given by Eg. 1121 In the 
top row a modelled GC harbouring 20 millisecond pulsars (A'^pgj^) 
is presented, while in the bottom row a model with A'^psr = 100 
is shown. In both cases, the distance to the cluster is taken to be 
d = 9.6 kpc. 



the ones produce by the clusters with larger pulsar popula- 
tion. 



4.3 ICS component to the 7-ray emission of the 
synthetic globular cluster 

The flrst step in obtaining an ICS emission component to 
the 7-ray radiation from globular clusters is creation of a 
cumulative spectrum of electrons injected by the population 
of synthetic MSPs into the GC environment. The electron 
ejection spectrum (see Sect. 13.4]) for a given pulsar (P, Bg, ce) 
is independent of a viewing angle, so once we have selected 
these basic pulsar parameters for all ApsR pulsars in the 
synthetic globular cluster, we construct the electron ejection 
spectrum for each synthetic pulsar: 



'I'cl,fe(7cl) = 
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where A/'d is the number of electrons with the energy _Eei = 
7oi TWcic^ ejected from pulsar magnetosphere per unit time. 
Then to obtain the cumulative spectrum of electrons, the 
electron ejection spectra from all the pulsars in the modelled 
globular cluster are summed up: 



^c 



= Yl *<=i.fe(7oi) 



(12) 



Examples of the cumulative electron spectra obtained 
for the population of millisecond pulsars residing in the syn- 
thetic globular cluster are presented in the right column of 
Fig. O The results for the cluster harbouring 20 synthetic 
millisecond pulsars is presented in the top right panel of 
Fig. [31 In the plot there are 25 different cumulative elec- 
tron spectra $oi displayed. A common feature present in 
all of the spectra is a narrow peak positioned at high elec- 
tron Lorentz factors. Its centre is around 7ei — 10®'^. A low 
energy tail is also characteristic of all the presented ^ci spec- 
tra. It stretches down to electron Lorentz factors ~ 10''. On 
average the particle flux level in the tail is two orders of 
magnitude lower than in the peak ("I>ei ~ 10^^ s~' MeV~' 
in the tail versus >I>ei ~ 10^^ s~' MeV"' in the peak). The 
shape of the tail is different between different globular clus- 
ter simulations. There are spectra where the tail is charac- 
terised by a type of plateau going from ■jei ~ 10^ to 10^, 
which is followed by a sharp decline towards lower Lorentz 
factor values. There are also cases where the low energy tail 
declines steadily towards small 701 values from the point it 
emerges from the high energy peak. The cumulative elec- 
tron spectra "l?ei obtained for the synthetic globular cluster 
harbouring 100 millisecond pulsars are presented in the bot- 
tom right panel of Fig. O There are 10 different spectra $ci 
presented in the plot. Similarly as for the case of the clus- 
ter harbouring 20 pulsars, the narrow high energy peak is 
present in all of the simulated electron spectra $ei- It is cen- 
tred at 7ci ~ 10^ ' 
flux is ~ 10^" ' 
low energy tail is by two orders of magnitude lower than 
in the peak. The shapes of the low energy tail are less di- 
verse with respect to the case of the synthetic cluster with 
smaller pulsar population. In all simulated cases the tail de- 
clines steadily towards low 7ci values right from the point it 
emerges from the high energy peak. The particle flux level 
in the tail changes from ~ 10^^ s~' MeV~' at its high en- 
ergy end (7ci ~ lO'^'*) down to ~ 10^^ s~' MeV"^ at its low 
energy end (7ei between ~ 10* and ~ 10^). 

In order to simulate the inverse Compton scattering 
component to the 7-ray spectrum of the synthetic globu- 
lar cluster, we are using a numerical model of the cluster 
developed by BS07. Because millisecond pulsars are con- 
centrated in the cluster core, it is assumed in the model 
that the relativistic electrons (their distribution is given by 
Eq. I12p are injected in the centre of the globular cluster, 
and further they diffuse gradually towards outer parts of 
the cluster. The diffusion of particles is treated in the Bohm 
diffusion approximation. On their outward way, electrons 
interact with the cluster magnetic field, and also with ambi- 
ent photon fields. BS07 (see fig. 1) showed that the radiation 
field within the GC and its nearby surrounding seems to be 
quite homogeneous. At the same time the electron diffusion 
distance mildly increases when moving towards the GC out- 



and the maximum value of the particle 
s^' MeV^'. The average fiux level in the 



skirts. Simple calculations show that the average diffusion 
distance for electrons injected in the cluster core and at its 
edge differ only by a factor of 2. Therefore the effect of dis- 
tribution of pulsars within the core of globular cluster on 
the cluster's ICS emission is expected not to be very signifi- 
canu (taking into account the uncertainties of the discussed 
scenario). 

Previously in their calculations, BS07 did not take into 
account synchrotron energy losses suffered by particles prop- 
agating through the cluster. In the presented calculations 
synchrotron losses are computed accordingly. However, the 
produced synchrotron radiation does not contribute to over- 
all 7-ray emission of the cluster. The characteristic energies 
of the synchrotron photons produced by electrons escaping 
from the pulsar magnetospheres and propagating within the 
globular cluster are far below the 7-ray range. They can be 
estimated from Schar ~ mdC^ (j?Gc/- Bcrit)7pi. where -Bcrit is 
the c ritical magnetic field (see e.g., iThompson fc DuncanI 
[1995I). For the electrons with Lorentz factors of the order of 
3 X lO'' and the magnetic field strengths within GC of the 
order of 10 ^G, the synchrotron photons have characteristic 
energies in the optical range, i.e. they will be very difficult 
to observe due to the huge background from the GC stars. 

In the calculations a following set of the globular cluster 
magnetic fields Bqc = (0.3, 1.0, 3.0, 10, 30) fj,G is taken into 
account. A main source of photons permeating the cluster 
are stars residing in the cluster itself. The stellar photons 
can be up-scattered by the relativistic electrons (with en- 
ergies up to ~ 1 TeV) via inverse Compton scattering up 
to TeV energies. The density of the stellar radiation field 
is not constant throughout the cluster, but it should mimic 
the distribution of stars in the system. For this reason, the 
density of stars in the clus ter is describ ed in the framework 
of the Michie-King model (|Michielll963l ): this allows for de- 
termination of the energy density of stellar photons t/gteiiar 
as a function of the distance from the cluster centre. In or- 
der to determine the energy density of the stellar photons 
within the cluster information on the cluster visual lumi- 
nosity I/vis is necessary. Moreover, the core radius Re, the 
half-mass radius Rh and the tidal radius Rt are also neces- 
sa ry for determining ?7ateiiar in the cluster (see eqs. 4 and 5 
of lBednarek fc Sitarekll2007h . 

As shown by BS07, the energy density of the stellar ra- 
diation field (t/stciiar ~ 300 eV cm^'^ inside a cluster core for 
the case where core radius is 0.5 pc and cluster luminosity is 
at the level of 10^ Lq) dominates over the energy density of 
the cosmic microwave background (CMB; Ucmb = 0.25 eV 
cm~'^). However, the latter photon field is important for esti- 
mating the ICS component (from tens MeV up to hundreds 
GeV) produced on the high energy electrons (energies above 
~ 1 TeV). The described processes are simulated using the 
Monte Carlo method, and lead to the production of the ICS 



•^ From the selected number of millisecond pulsars used 
to model synthetic 7-ray spectra of 47 Tuc and Ter 
5 (see Sect. |5) only 5 out of 23 and 1 out of 31 
l |http://www.naic.edu/~pfreire/GCpsr.html| , respectively, lie 
beyond the distance of 2 x _Rc but within the radius of 3 x i?c (^c 
is the cluster core radius). This is relatively small number. Thus, 
assuming that all MSPs are situated in the cluster core does not 
cause the ICS emission to be overestimated. 
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emission spectrum of the globular cluster. The details of the 
numerical code are presented in the work of BS07. 

The results of the numerical simulations performed for 
two globular clusters: Terzan 5 and 47 Tucanae are presented 
in Sect.O The outcome of the simulations is the total 7-ray 
spectrum for each of the clusters. The spectrum consists of 
the magnetospheric and ICS component. 



5 COMPARISON OF THE SIMULATIONS 
WITH THE OBSERVATIONS 

Using the numerical procedure presented in Sect. [4l 7-ray 
spectra were computed for two selected globular clusters: 
Terzan 5 and 47 Tucanae. Where possible, the simulated 
spectra are compared with the publicly available observa- 
tional data. 

An important comment about the selection of pulsar pe- 
riods for the population of MSPs residing in the clusters is 
necessary. Because for each of the selected clu sters a number 
of millisecon d pulsars w as detected (sec e.g., Ranso m et al.l 
^005; Manc hester et all , 1991: Caniilo ct al. 2000) and their 
spin periods were determined, a different approach to the 
problem of the spin period P selection than the one pre- 
sented in Sect. [4] was used. From the list of the known MSPs 
in each cluster, pulsars with periods satisfying the criteria 
P < 11 ms were selected as the sole members of the popula- 
tion of the millisecond pulsars in the cluster. Knowing their 
spin periods and the density of the simulated grid in our 
database (Sect.[33Ji for each MSP the spin period value as 
close as possible to the true P value was assigned from the 
database. For each MSP the value of the magnetic field Bs 
and the inclination angle a was randomly selected from the 
discrete values simulated in our database (Sect. [373)) . From 
this point onwards, the procedure leading to the production 
of the cluster 7-ray spectrum matches the one presented in 
detail in Sect.]!] 



5.1 The case of Terzan 5 

Currently, Terzan 5 is the globular cluster with the largest 
number of the detected millisecond pulsars. There are 34 
MSPs observed in the cluster (see the catalogue of P. C. 
Freirqj). In addition, it is one of the f irst globular clus- 
ters detected in 7-rays by Fermi/LAT (|Abdo et al.l I2OI0I : 
iKong et al.ll20ld ). Its luminosity estimated in the energy 
range between 100 MeV up to around 20 GeV is L^ ~ 
2.6 X 10^^ erg s"^ (|Abdo et all I2OIOI ) assuming the clus- 
ter is situated at a distance of d ~ 5.5 kpc from the Sun. 
The observed 7-ray spectrum can be well fitted with a 
power law with an exponential cut-off N{£) ~ S^^e^ ' " 
{N{£) = dN-y/d£ is a photon flux, F is a photon index, and 
£c is a cut-off energy). The spectral pa rameters estimate d 
for Ter 5 are F ~ 1.4 and fc ~ 2.6 GeV (|Abdo et al.ll2010D . 
Recen tly, the cluster has been d etected by the H.E.S.S. tele- 
scope (|Abramowski et al.ll201ll ) in the TeV domain. The ob- 
served photon flux in the energy range above 0.4 TeV is 
~ 1.2 X 10~^^ cm~^ s~^, while the photon index of the spec- 
trum fitted with the power law N{£) ~ £^^ is F ~ 2.5. The 
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TeV emission detected from the direction of Ter 5 is slightly 
shifted with respect to the clust er core and it extends wel l 
beyond the cluster tidal radius. lAbramowski et al.l (|201ll ) 
speculate that the observed TeV emission may be a source 
coincidence, e.g. with a PWN associated with a radio-quiet 
pulsar, or that the TeV 7-rays could have hadronic origin. 
If neither of these scenarios is true, then Ter 5 is the only 
globular cluster detected so far in the TeV domain with the 
Cherenkov telescopes. 

Important parameters of Ter 5 cluster that are used in 
calculations of the ICS spectral component are the visua l 
luminosity of the cluster Lvis — 1-5 x 10^ Lq (|Harrij|l996l ). 
The core radius Re, the h alf-mas s radi us _Rh and the tidal 
radius _Rt were taken from iHarrij (|l996l ). The ambient pho- 
ton fields up scattered via inverse Compton process are the 
cosmic microwave background (C/cmb = 0.25 eV cm~'^) and 
the photons from the stellar population in the cluster (the 
energy density was estimated with eqs. 4 and 5 of BS07). 

Out of 34 known MSPs in the cluster, 3 have spin pe- 
riods larger than 11 ms. Thus, A^psr = 31 were used to 
simulate 7-ray radiation from the cluster. The results are 
presented in Fig. 2] For one set of NpsB., with chosen (P, Bs, 
a) one cumulative spectrum of electrons (given by Eq. I12|) 
is obtained. This spectrum is then injected into the clus- 
ter environment and it interacts with the ambient photon 
fields leading to the production of the ICS component of 
the 7-ray spectrum of the modelled globular cluster. Simul- 
taneously, for one set of A^psr many different cumulative 
magnetospheric spectra (given by Eg llOl) are simulated. The 
cumulative magnetospheric spectrum is the component of 
the synthetic cluster 7-ray spectrum stretching from around 
MeV up to around few tens of GeV. The grey area in the 
bottom panel of Fig. [4] shows the region occupied by the 
group of different magnetospheric spectral components sim- 
ulated for the cluster. It can be treated as the uncertainty 
of the GC emission spectrum in the MeV to GeV energy 
domain. In the bottom panel of Fig. U the simulated total 
7-ray spectra for Ter 5 are presented with solid lines. To 
construct the total spectra, one magnetospheric component 
out of the simulated group was chosen. Five different ICS 
components were obtained for the cluster. They were sim- 
ulated taking different values of the cluster magnetic field 
Bqc (see Sect. 23]). For clarity, these ICS components are 
presented in the top panel of Fig. [4] separated from the cho- 
sen magnetospheric component. 

The overall shape of the ICS spectral component (see 
the top panel of Fig. [4]) can be described as having two peaks: 
the high energy peak (at photon energies above 1 TeV), and 
the low energy peak (at f ~ 10 GeV). In the presented to- 
tal GC spectra (the bottom, panel of Fig. |4| only the high 
energy peak is visible, the low energy peak is dominated by 
the fading magnetospheric component of the cluster 7-ray 
emission. The ICS spectral component has slightly different 
shape and also different level depending on the cluster mag- 
netic field values Bgc used in calculations (see Sect. 14. 3p . 
For low magnetic field strength in the cluster (Bgc ^ 1 mG) 
a clear dip between the high energy peak (resulting from 
IC scattering of the starlight photon field) and the low en- 
ergy peak (resulting from IC scattering of the CMB) is vis- 
ible. The dip gradually disappears when the magnetic filed 
in the cluster increases (-Bgc > 1 A^G), leading to single 
peaked ICS spectrum for the case where Bgc = 10 and 
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through the globular cluster environment can be explained 
by the synchrotron loses suffered by the electrons in the 
high magnetic fields. For low values of the magnetic field 
_Bgc in the cluster, diffusion of electrons occurs rather fast 
(fdiff — 3 X 10^^ s for particles with energies E^i = 0.1 TeV, 
the cluster magnetic field Bgc = 1 mG and the cluster half 
mass radius 7?h = 1 pc) and their IC cooling is less effi- 
cient. For strong magnetic fields particle diffusion is much 
slower (idiff — 9 x 10^^ s for Bgc = 30 /^G). However, as the 
-Bgc increases, the synchrotron loses of the particles become 
larger leading to the saturation and decrease in production 
of the ICS photons. 

In Fig. 2] the results of simulation for Ter 5 are also 
contrasted with the available observational da ta in the 
MeV to GeV (F ermi/LAT data:lAbdo et al.ll2010l ) and TeV 
(H.E.S.S. data; lAbramowski et al.l |201l|) energy domain. 
With our simulations we obtain the magnetospheric 7-ray 
spectra that may reproduce or underestimate the cluster 
emission observed at energies below 10 GeV (see the results 
presented in the bottom panel of Fig. |4]). However, we are 
unable to reproduce the exact spectral shape observed at 
TeV energies. 

The 7-ray flux calculated in this paper should be consid- 
ered as a guaranteed lower limit provided that the model for 
high energy processes occurring in the inner pulsar magneto- 
sphere is correct. We conclude that TeV 7-ray observations 
of GCs can provide independent constraints on the pulsar 
models. 



Figure 4. Synthetic spectra calculated for the globular cluster 
Terzan 5. The total number of millisecond pulsars in the clus- 
ter is A^psR = 31. The distance to the cluster is assumed to be 
d = 5.5 kpc. Top panel: ICS spectra computed for different clus- 
ter magnetic field values Bqq. Solid black, green dotted, green 
dashed, green dot-dashed and brown dot-dot-dashed lines show 
models with Sgc of 0.3, 1.0, 3.0, 10, 30 /xG, respectively. Bottom 
panel: Total 7-ray spectra (including magnetospheric and ICS 
component) are depicted with solid lines. Colours of the lines 
correspond to the value of the magnetic field in the globular clus- 
ter - from black, through shades of green to brown correspond to 
values of -Bgc from 0.3 fj,G up to 30 fiG. The colour-coding of the 
lines is the same as for the top panel. ICS spectra in the bottom 
panel were re-binned to match the resolution of the magneto- 
spheric component spectra. Grey area shows the region occupied 
by the magnetospheric spectra of the cluster obtained for pre- 
sented simulation runs (this can be treated as the uncertainty 
of the MeV to GeV emission spectrum of GC). Triangles show 
the 7-ray spectrum (with u ncertainties) of Te r 5 as observed by 
the Fermi/LAT instrument IIAbdo et al.ll2010l) . Diamonds depict 
spectrum (and its errors) obser ved from the direction of the clus- 
ter by the H.E.S.S. telescope llAbramowski et 311120111 ). Dashed 
line shows the CTA differen tial sensitivity curve with 50 h of in- 
tegration llActis et al.ll201ir) . 



30 /iG. When we increase the -Bgc value in the cluster, we 
see that the level of the ICS emission of the cluster increases. 
For the magnetic field values around 3 /iG the emission 
reaches maximum. Further increase in Bgc results in the 
decrease of the emission level of the ICS component. This ef- 
fect of the increase-saturation-decrease of the emission com- 
ponent resulting from the up-scattering of the stellar pho- 
ton field and the CMB on the relativistic electrons diffusing 



5.2 The case of 47 Tucanae 

The 47 Tucanae globular cluster harbours 23 millisecond 
pulsars de tected so far through radio searc hes (see e.g., 
ICamilo fc Rasig 2005: Manchest er et al.lll99ll '). Similarly to 
Ter 5, 47 Tuc has been detected in the MeV-GeV energy 



domain with the Fermi/LAT instrument (JAbdo et al.ll201C 
Its 7-ray spectrum can be approximated with the exponen- 
tially cut-off power law with a spectral index F ~ 1.4 and a 
cut-off energy £c — 2.2 GeV. The cluster luminosity in this 
energy range is L^, ~ 4.8 x 10^* erg s~^ for the cluster dis- 
tance d of 4.0 kpc. No TeV emission has been reported so far 
from 47 Tuc. However, the u pper limit to the TeV e mission 
of the cluster was estimated (jAharonian et al.ll2009l ). 

Th e 47 Tuc clu ster visual luminosity is Lvis — 7.5 x 
10^1/0 (|Harrislll996l ). The core radius Re, the half-m ass ra- 
dius i ?h and the tidal radius Rt were taken from iHarrisI 
([199a). The ambient photon fields up scattered via inverse 
Compton process are the cosmic microwave background 
(C^cMB = 0.25 eV cm~^) and the photons from the stellar 
population in the cluster (the energy density was estimated 
with eqs. 4 and 5 of BS07). 

All of the detected 23 MSPs in 47 Tuc have spin periods 
smaller than 11 ms, so A^psr = 23 in the simulations. The 
results are presented in Fig. [5] Similarly as for Ter 5, only 
high energy peak of the ICS component is visible in the clus- 
ter total spectrum (see the bottom panel of Fig.[S}. The low 
energy peak, which is visible in the ICS spectra presented in 
the top panel of Fig. O is dominated by the magnetospheric 
component. The increase-saturation-decrease behaviour of 
the ICS emission component is visible, though the level dif- 
ference between the ICS spectrum obtained for the lowest 
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Figure 5. Same as Fig. |4] but for tlie globular cluster 47 Tu- 
canae. Blue dot-dashed line represents the upper limit on TeV 
emissi on of 47 Tuc as obtain ed from H.E.S.S. telescope observa- 
tions iJAharonian et al.ll200!J ). Black dashed line sho ws the CTA 
differ ential sensitivity curve with 50 h of integration IIActis et alj 

MB). 



6 CONCLUSIONS 

The 7-ray spectra (energies above ~ 100 MeV up to ~ 
10 TeV) for synthetic globular clusters were calculated. 
The magnetospheric component (photon energies between 
~ 100 MeV and ~ 30 GeV) of the cluster emission results 
from the cumulative 7-ray emission of millisecond pulsars 
residing within the cluster core. The ICS component (pho- 
ton energies above ~ 30 GeV) is produced by the ambi- 
ent photons (the stellar photon field originating from stars 
within the cluster, and the cosmic microwave background) 
up-scattered via inverse Compton process on relativistic 
electrons injected into the cluster environment from the syn- 
thetic MSP magnetospheres. The detailed simulations were 
performed for two globular clusters, Terzan 5 and 47 Tu- 
canae, for which observational data are available in the en- 
ergy range between MeV and TeV. Our calculations are able 
to reproduce, in both cases, the spectrum obser ved with the 
Fermi/LAT instrument for energies < 10 GeV (JAbdo et al-l 
I2009al . l2010l ). However, for Ter 5 where TeV observations are 
available, our simulations cannot account for the observed 
spectral shape in this energy domain. 

On the other hand, it is possible that colliding pulsar 
winds, as proposed by BS07, can accelerate leptons to rel- 
ativistic energies. Note that other com pact objects wit hin 
the GCs, such as rotating white dwarfs (|Bednarekll2012l ) or 
accreting neutron stars, could be responsible for accelera- 
tion of electrons and their injection into the volume of GCs. 
The energy spectra of these shock-accelerated electrons are 
power-law in nature. Inclusion of this additional source of 
relativistic leptons - of different spectral character than the 
ones produced in MSP magnetospheres - into the calcula- 
tions of ICS component could improve the match between 
the synthetic TeV spe ctrum calculated for Te r 5 and the 
H.E.S.S. observations (JAbramowski et al.ll201ll ). Thus, the 
calculated synthetic spectra of globular clusters at TeV en- 
ergies should be treated as a lower limit to GCs emission 
in this energy range. This lower limit can be fu rther tested 
by th e future Cherenkov telescopes like CTA (JActis et al.l 
l20TIh . 



value of the magnetic field Bgc = 0.3 /iG and the spectra 
obtained for higher Bgc is smaller than in the case of Ter 5. 
The simulated 7-ray spectra are compared with the 
availa ble observational data from the Fermi/LAT instru- 
ment (jAbdo et al.ll20ld V Our results rather nicely reproduce 
the cluster spectrum observed for energies < 10 GeV. As the 
cluster has not been det ected in the TeV dom ain, we plot 
the H.E.S.S. upper limit (|Aharonian et al.ll2009l '). The mod- 
elled ICS spectral components do not violate the H.E.S.S. 
upper limit. The modelled TeV emission of 47 Tuc falls right 
below the presented upper limit. The differentia l sensitivity 
curve for the Cherenkov Telescope Array (CTA: lActis et al.l 
I2OIII ). the next generation system of Cherenkov telescopes, 
is also displayed in Fig. [5] The curve shows the CTA sen- 
sitivity obtained after 50 h of integration. Our simulations 
show that the cluster should be detectable with the CTA 
for all of the studied Bgc- However, the ICS emission is go- 
ing to be more pronounced for the moderate values of the 
magnetic field within the cluster (Bgc between 1 fiG and 
10 ^G). 
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